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Abstract 

In this work we investigate aspects of light cones in a Schwarzschild 
geometry, making connections to gravitational lensing theory and to a 
new approach to General Relativity, the Null Surface Formulation. By 
integrating the null geodesics of our model, we obtain the light cone 
from every space-time point. We study three applications of the light 
cones. First, by taking the intersection of the light cone from each 
point in the space-time with null infinity, we obtain the light cone 
cut function, a four parameter family of cuts of null infinity, which is 
central to the Null Surface Formulation. We examine the singularity 
structure of the cut function. Second, we give the exact gravitational 
lens equations, and their specialization to the Einstein Ring. Third, as 
an application of the cut function, we show that the recently introduced 
coordinate System, the "pseudo-Minkowski" coordinates, are a valid 
coordinate System for the space-time. 



1 Introduction 

The purpose of this work is to develop a simple model space-time in which we 
study gravitational lensing and the light cone cuts of null infinity. These cuts 
are central to a recent reformulation of General Relativity known as the Null 
Surface Formulation The model we consider consists of a Schwarzschild 

exterior region surrounding a spherically Symmetrie, static, constant density 
dust star. 
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The Null Surface Formulation makes explicit use of future null infinity, 
denoted by X + , which has the topology of R x S 2 . In an asymptotically 
simple space-time, X + represents the future endpoints of all null geodesics 
and can be added as a boundary to the physical space-time through a process 
of conformal compactification ||. The Standard coordinates on X + are the 
Bondi coordinates (u, 9, 0), where u labels the R part, and 9 and <fi label the 
sphere. 

On X + , a light cone cut is the intersection of the light cone from a partic- 
ular point in the inferior with future null infinity. In Bondi coordinates on 
X + , the light cone cuts are given by a cut function, 

u = Z(x a ,9,(f)), (1) 

where for a fixed initial point, x", the cut is a deformed sphere, possibly with 
self-intersections and singularities. The crux of the Null Surface Formulation 
is that from the knowledge of the light cone cuts, Eq. (HD, one can reconstruct 
all of the conformal information of the space-time. 

The current paper explores the light cones from an arbitrary point in 
a Schwarzschild space-time surrounding a constant density, dust star. By 
integrating the null geodesics in an inverse radial coordinate l = l/v2r, we 
obtain parametric expressions for the future and past light cones of an initial 
point, x", in terms of a "distance" /, and two "directional" parameters which 
span the sphere of null directions at x®. The intersection of the future light 
cone of an arbitrary initial point with null infinity gives the light cone cuts. 
In Section 3, we study the singularity structure of the cuts, explaining how 
the singularities are related to the formation of conjugate points on the light 
cone. We show in Section 4 that the equations for the past light cone obtained 
in Section 2 are, in fact, examples of exact lens equations. As a special case, 
we give the exact formula for the Observation angle for an Einstein Ring. In 
the final section, we show that the "pseudo- Minkowski" coordinates, defined 
in ||, form a valid coordinate System for the entire space-time. 

2 Integrating the Null Geodesics 

To begin, we integrate the null geodesics of a Schwarzschild space-time with 
an inferior constant density matter region. Since we eventually consider each 
geodesic's limiting endpoint at X + in order to obtain a cut function, we in- 
tegrate the null geodesics using a conformal Schwarzschild metric which is 
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regulär at null infinity. The Integration is performed using a "radial" Param- 
eter l = l/(\/2r), so that l = corresponds to the point at null infinity, 
while a finite l > will be a point in the interior. 

The general form for a static, spherically Symmetrie metric with signature 
(+, -, -, -) is given by 

ds 2 = f(r)dt 2 - h(r)dr 2 - rW, (2) 

where r 2 dQ 2 = r 2 d9 2 + r 2 sin 2 Odcf 2 is the line element on the sphere. In 
our model, the funetions f{r) and h(r) will be continuous, pieeewise smooth 
funetions for an exterior Schwarzschild region and an interior constant density 
dust Solution with a radius R: 



f(r) = 1 - — = f ext , r > R, 



h ( r ) = ~, 2M^ = ^nt, r < R, 

1 R 3 

K r ) = 1 ZTW. ^ hexU r > R - (3) 

r 

We assume that the radius of the interior region extends beyond the r = 3M 
unstable circular orbit for null geodesics, ensuring that the space-time is 
asymptotically simple. Working in a retarded null coordinate, u = t — 
J dr J h{r) / f (r) , and the inverted radial coordinate, / = l/(\/2r), a confor- 
mally rescaled version of the metric, Eq. (|2|), which is regulär at null infinity, 
is 



ds 2 = AI 2 f(l) du 2 - A^Jk(l)du dl - dü 2 , (4) 
where r is replaced in terms of the variable l and k{l) is given by 

k(l) = f(l) h(l). 
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c = 21 2 f(i) ü 2 - 2^k[T) ü i - -e 2 , (5) 



It is convenient to integrate the null geodesics first in a plane, and then to 
use the spherical symmetry to rotate the Solution to an arbitrary orientation. 
The restriction we employ is to take a particular initial point, x%, lying on 
the —z axis, 

%o = ( u o, lo, o = tt, <p = 0), 

and the geodesics as lying in the x-z plane. To dehne the x-z plane, one 
usually allows 9 to ränge from to tt and has <fi = or tt. In order to 
facilitate our discussion, we will not use this Convention. Instead, we require 
that = and allow a variable, 9, to ränge from to 2n. In terms of the 
variables (u, l, 0), a Lagrangian corresponding to the conformal metric is 

ir 

2 

where dots indicate derivatives with respect to an affine parameter r. The 
geodesic equations are 

^(2/ 2 /(o « - V*(Ö i) = o, 

^(-V^Ö «) - (2/ f(l) ü 2 + l 2 /'(/) ü 2 ) + ^=W)) = 0, 

A = o, (6) 

with primes denoting derivatives with respect to /. To solve for null geodesics, 
we impose the null condition on the Lagrangian, 

£ = 2l 2 f(l) ü 2 - 2^k[T) ü 1 - -6 2 = 0, (7) 

and use this equation and the hrst and third geodesic equations as inde- 
pendent equations for (w,/,0). After hnding some trivial hrst integrals and 
rearranging, independent equations for null geodesics can be written as 



1 + Jk(l) l 

11 = 

2/2/(0 ; 

h - h 2 

i = ± 



i-b 2 p f(i) _ r— 
\ — W) — = ± v A(/) ' 
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e = b. (8) 

In these equations the sign of l indicates whether the geodesic is incoming 
(positive) or outgoing (negative), and the parameter b is a free integration 
Parameter labeling the initial direction of the geodesic. In Appendix A, 
we show how the parameter b is related to the initial direction of the null 
geodesic. 

In integrating the equations, geodesics which are initially outgoing will 
have b values ranging from zero, corresponding to a ray traveling radially 
outward, to some maximum value b m , for which the geodesic is initially tan- 
gent to a sphere of radius r Q = l/(\/2l ). The value of b m is the value of b 
for which / = at the point l = l - 

(9) 



For geodesics which are initially incoming, the ränge in b will be from b m back 
down to zero. The value of l will increase until reaching some maximum value, 
l = lb, which is a minimum r value. The value of lb is the Single (real) root 
of the equation, l = 0, or the root of: 

1 - b 2 l\ f(l b ) = 0. (10) 

At lb, l changes sign, and the geodesics will head out to infinity. 

It is convenient to reparametrize the equations of motion, Eqs. (||), us- 
ing l instead of the affine parameter, and to express the Solution to the null 
geodesic equations in terms of integrals over /. In terms of these integrals, 
geodesics on the light cone connecting the initial point, x" = (w , / D , 9 a = ir, (p Q = 0), 
with the final point, x a = (u, l, 6, = 0), are given by 



u = u(u Q , Z , l, b) = u Q + (— l) e / dl' 



1 ± y/kW A{V) 
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where the integrals are defined piecewise over the various segments in / and 
k, and the appropriate signs are chosen if the geodesics are incoming (+) or 
outgoing (— ). The future light cone is given by e = 0, and the past light 
cone by e = 1. 

The füll Solution to the null geodesic equations is obtained by performing 
a rigid rotation of the spatial plane defined by the path of this geodesic and 
the z axis, to an arbitrary orientation. Due to spherical symmetry, the light 
cone for an arbitrary point will be axially Symmetrie about the spatial line 
connecting that point and the spatial origin (defined by r = 0). We will call 
the angle of revolution about the axis of symmetry 7; the angle 7 essentially 
defines an orientation of the spatial plane containing the geodesic and the 
axis of symmetry. (In the case that the initial point lies on the —z axis, 7 
is simply the angle 0.) To obtain the füll light cone from an arbitrary initial 
point, we need to rotate the Solution that is Symmetrie about the z axis over 
to the axis of symmetry defined by the arbitrary point. We can think of this 
as rotating the particular initial point x a a = (u , l , O = n, O = 0) to 
the general initial point x% = (u , l , d , <f> ) and allowing the orientation 
Parameter 7 to take any value between and 2n. This rotation is explicitly 
performed in Appendix B. 

It is often convenient to use complex stereographic coordinates, defined 



instead of the angular coordinates 6 and 0. Throughout this paper we freely 
switch back and forth between the two coordinate Systems. An arbitrary 
initial point is given in stereographic coordinates by x" = (u , l a , ( Q , ( Q ). 

The füll light cone, obtained by using the rotation of the Solutions of 
the null geodesic equations in the x-z plane to an arbitrary initial point and 
orientation given in Appendix B, are expressed parametrically in terms of 
the parameter / as 



by 




(12) 



u = u(u , l , l, b) 




± 



1 ± y/k(V) A{V) 



) 



■o 



4Ww 2 /(*')) 



C(k,Co,Co, J,&,t) 



( 



( 



e f7 CQt e(^) + v /^- e - t7 



) 



6 



C(^O) Co 5 Co; b, 



7) 



( 



c 
c 



o 




e -f7 CQt e(^) + ^ ei7 



(13) 



o 



v /^ e -f7 cot e(^) +ei7 



where 0(1, l D , b) is given by the integral in Eq. (JTTJ) 




(14) 



and the Convention for e and signs are taken as before. 

Eqs. ( |i~3| ) represent the entire future and past light cones for an arbitrary 
point in the space-time. The light cones are given in terms of two parameters, 
(b, 7), which span the sphere of initial null directions at the initial point 
x". We note that the u coordinate does not depend on 7 due to the axial 
symmetry of the light cone. 

3 Light Cone Cuts and their Singularities 

While the light cone from an arbitrary point in Minkowski space is always 
smooth, light cones in an asymptotically simple space-time have, in general, 
self-intersections and several different kinds of singularities. These singu- 
larities are directly related to the formation of conjugate points along null 
geodesics. A pictorial representation of a light cone with singularities is given 
in Figure [I]. Since the light cone cut function is the intersection of the fu- 
ture light cone with null infinity, it inherits the singularity structure of the 
light cone. In this section, we study the cut function in our model as a 
representative example of the Null Surface Formulation. 

A parametric version of the cut function is obtained by setting the values 
of / and e to zero in Eqs. flTB]) : 




(15) 




(16) 
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uuMD-m- 4^±ä1|, (17) 

where the function Q(l ,b) is given by 



e(/ ,6) = 7r -/° ^ 



(18) 



If it were possible to invert the pair of equations, Eqs. (|iq) and (|17D , for 6 
and 7, obtaining the functions 

"7 ^-*(Coo, Coo, Co; Co; ^o), (^) 

and 

6 = ß(Coo,Coo,Co,Co,/o), (20) 

one could produce the füll cut function in the form of Eq. ([!]), 

u = Z(x a ,(,(), 

by inserting the Solution for b in terms of (£, £) from Eq. ( f2Ö| ) into the Mqo 
Solution, Eq. (1151) . 

While Coo(/o, Co, Co, b, 7), and Coo(^o, Co, Co, b, 7) are Single valued 
in b, they will not, in general, have unique inverses — more than one initial 
direction acquires the same value of C or C at 1+. This implies that there 
are singularities in the cut function itself and that the global inversion of 
Eqs. (|i~6|) and (|17|) for b and 7 will be impossible. In such a case, we will not 



be able to find an explicit cut function in the form of Eq. ([[]), but will be 
forced to work with the cut function in a parametric form. 

In a sense, singularities in the light cone cuts are places where the Null 
Surface Formulation undergoes technical difficulties — a natural coordinate 
System used in the theory is not well defined at these points. We now be- 
lieve that these difficulties can be overcome by using a particular parametric 
representation the cut function. A primary interest here is to study the 
singularities in the cut function. 

There is a complete Classification of the stable singularities of the cut 
function which can be applied to our model, due to Arnol'd and his collab- 
orators || ||. A stable singularity is one which does not disappear under 
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small perturbations. For two dimensional surfaces, such as the light cone 
cuts obtained by fixing the initial point in the cut function, there are only 
two types of stable singularities. These are the cusp ridge and the swallow- 
tail. Due to high level of symmetry in our model, any cut function must 
be axially Symmetrie. This means that, although the cut function is a two 
dimensional surface, it can be represented by a one dimensional curve whose 
revolution about some axis gives the cut function. For a one dimensional 
curve, the only stable singularity is a cusp, which implies that we will not 
see swallowtail singularities in our model. 

In Figure ||, we give plots of three cut funetions, suppressing the axially 
Symmetrie dimension, for three different values of the initial radial Posi- 
tion, / . These figures show that singularities appear as the initial point 
moves away from the spatial center of the space-time. A smooth cut of null 
infinity, corresponding to a cut from the light cone of a point close to the 
center of the space-time, will be a smooth sphere-like surface. Due to the 
axial symmetry, a cut with singularities will have a circular cusp ridge and a 
single crossover point. Figure [| gives a pictorial representations of a singular 
cut. 

Because the cusp ridge singularity in our cut function is stable, it repre- 
sents a generic possibility for a cut function in a general space-time. Specif- 
ically, this singularity will remain if one makes a small perturbation of the 
metric away from a Schwarzschild metric. The crossover point along the cut 
function is an unstable singularity, arising from the high degree of symmetry 
in the Schwarzschild case, and is directly related to Einstein's Rings, the 
astronomical phenomenon where a spherical lens causes lensing in a uniform 
circular shape |7|]. We discuss the rings in the next section on gravitational 
lensing. 

The singularities in the cut function represent points which are conjugate 
to the initial point x a . For a fixed x", we can think of the parametric equa- 
tions for the cut function, Eqs. (p~5|) , fll6|), and (|17|), as a map between the 
initial directions of the geodesics at the initial point and the final position 
(«oo, Coo, Coo) at T + . One way to find the points of X + which are conjugate to 
the initial point is to find the points at X + for which the Jacobian matrix 
expressing the mapping, 
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J 



V 



db 

dCoo 
db 

db 



dC,oo 
dj 

ÖCoo 
97 



(21) 



/ 



drops in rank |J. For this to occur, all three 2x2 determinants must be 
zero. 

With no loss in generality, we consider the cut function of an initial point 
on the —z axis. In this case, the cut function is obtained by setting the final 
value of l in Eq. (|TT|) to zero, and restoring the rotational degree of freedom 
by setting 7 = 0oo. Thus, the cut function for an initial point on the —z axis 
is written in terms of the coordinates («oo, Ooo, 0oo) as 



Woo(Wo, lo, b) = U Q + 



dl' 



1 ± Jk{l>) A(l>) 



M(Z')(2Z' 2 /(//)) 



0<x>(la,b) = 71 



dl' 



A{V) 



0oo = 7- (22) 
In this case, the Jacobian matrix corresponding to Eq. (pl| ) is given by 

/%* 0' 



J 







db 

V 1 



(23) 



It is clear that for the Jacobian to drop rank we must have 



and 



döoo 
q = = 0, 



db 



(24) 



(25) 



A combination of numerical and analytic calculations shows that these two 
conditions are satisfied simultaneously only at the cusp ridge shown in the 
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figures, and hence the cusps on the cut function are conjugate points to the 
initial point. 

An alternative way of deducing the singular points is to consider the first 
and second derivatives of Uoo with respect to 9oo- For a cusp singularity, 
the first derivatives are always finite, while the second derivatives diverge. 
Working parametrically in b, these derivatives are 



du r 



d6 n 



db 

Ö9oc 

db 



and, 



d 2 u no ( 86, 



-2 r 



oo 



oei, v db 



d 2 u^ du„o d 2 e c 



db 2 de^ db 2 



(26) 



Numerical computation shows that the first derivative is finite and non-zero 
along the cusps in the light cone cut. As written, the second derivative 
is indeterminant along the cusps since the term in Square brackets is zero 
there. Applying L'Hospital's rule, one sees that the second derivative actually 
diverges as q^ 1 . 

As a final note on the cusp singularities, we list the behavior of sev- 
eral important quantities in the Null Surface Formulation as one approaches 
the cusps along the cut function. All of these quantities are computed as 
derivatives of the cut function with respect to the complex stereographic 
coordinates ( and (, and have been computed in this model parametrically 
for initial points along the — z axis. These derivatives are denoted by 5 and 
9 §. In terms of q, which approaches zero as one approaches the cusp, 
the behavior of some important quantities of interest are listed below for 

Quantity Behavior 



üü = bZ regulär 
A = Ö 2 Z q- 1 



r = mz q- 1 

|A,i| = l§l 1 
Q 2 = g ab Z, a mz 
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4 Gravitational Lensing Equations 



An important goal of gravitational lensing theory is to construct lens equa- 
tions which give the position of sources in terms of directions seen by an 
observer and distances to the source. Typically, lens equations are obtained 
via approximations on the kinematics of the null geodesics of the source 0. 

Recently, a way to produce completely general, exact lensing equations 
has been found, and a paper is being prepared which develops gravitational 
lensing theory from this perspective |J. Our model provides an explicit ex- 
ample of such a formulation. In this section we give exact lensing equations 
for the Schwarzschild space-time with a constant density dust inferior re- 
gion, and show that our exact equations reduce to Standard approximate 
lens equations. In the special case that the source, lens, and observer are 
spatially colinear, we give an exact expression for the Observation angle in 
an Einstein Ring. 

In Section 2, we derived the future and past light cones of an arbitrary 
point in the space-time. Recall that a lens equation should express the lo- 
cation of the source in terms of some "distance" from the observer, and the 
directions which the observer views the geodesics on the past light cone. The 
equations of the past light cone, Eqs. (|[|), are such a set of equations. In 



these equations, the observed directions are given by the parameters (6,7), 
and / gives the "distance." Hence, exact lens equations for our model are 



C(/ ,Co,Co,/,&*,7* 




e i 7 * CQt e^T) + y^e-K 



VW V-VCoCoeT-r* cot + e *' 



with 



e(i,i ,b*) =tt- [ l di' \ ±-= 



-b* 2 V 2 f(l') i 



In these lens equations, the spatial location of the source is the point 
(l, C,C)- For an observer at the point (7 ,Co,Co), the observed directions of 
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the geodesic on the past null cone are given by the particular values of (b, 7) 
which connect the source and the observer, denoted as (6*, 7*). There may, 
in fact, be more than one set of values for (6*, 7*), as the process of focusing 
may produce more than one "image." 

Due to spherical symmetry, any observer may be considered as lying on 
the —z axis, and the source may be taken as lying in the x-z plane. We are 
interested in the case where the lens is situated between the observer and 
the source, and when the rays do not pass through the inferior region of the 
star. In this case, the lens equations, Eqs. (|27|), reduce to the Single equation 
for 6 which was found in Eq. (O): 



0(/ o , l,b*)= TT- / dl 1 ±-== (28) 

-b* 2 l' 2 /(/') 



This lens equation specifies the location of the source in terms of the observed 
direction of the geodesic, b*, and a "distance," /, to the source. 

In Appendix A, we find the relationship between the angle at which a null 
geodesic crosses the z axis, the parameter b, and a position l. In the lensing 
case, this "Observation angle," denoted by ipobs, is related to the observer 
position, l a , and the observed direction, b*, by 

,* sin^ (29) 



By replacing b* by ip b s in the lens equation, Eq. (|28]) , the lens equation 
takes a more conventional form, where the direction parameter is the actual 
Observation angle: 



Q(l Q ,l,b*)=it 



lo Jf(lo) 



dl' 



1 



V 



, _ sm 2 j> obs 1'2 f(l>) 

L WJW) 



(30) 



A typical approximate lens equation for the Schwarzschild model [[<]] is 

2R S D LS 



(31) 



In this approximation, ß is the Euclidean angle between the source and the 
center of the space-time, and Rs = 2M is the Schwarzschild radius. The 
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Euclidean distances between the source and lens, source and observer, and 
lens and observer, are given by -Dls, Ds, and Dl respectively. Figure f| shows 
the case under consideration. We now show that our lens equation, Eq. fl3ÖD 
or Eq. (|2~8D, reduces to the approximate formula, Eq. (|31~D, under appropriate 
approximations. 

Taking into account the correct signs for incoming and outgoing rays, the 
right hand side of Eq. (|28|) can be written as 

= TT- A(M,b*,k,l), (32) 

where 



b* dl' 



+ 



2V2M b* 2 l' 3 - b* 2 V 2 + 1 
b* dl' 

2V2M b* 2 l' 3 - b* 2 l' 2 + 1 



(33) 



Here, / D is the position of the observer, l is the position of the source, and 
lb is the value of / for which the geodesic comes closest to the lens, attained 
when / = 0. The maximum l value, £&, is, from Eq. (|IÖD , the Solution of the 
equation 

1 - b* 2 Ii + 2V2M b* 2 l 3 b = 0. (34) 

For convenience, we assume that the source is closer to the lens than 
the observer, so that l > l - To proceed, we assume that the dimensionless 
quantities M l = e and M l Q < e are small and make a Taylor series expansion 
of A in terms of e: 



A(e,&V o ,0 = A(e = 0,&V o ,Z)+e 



dA 
~de 



= A + Ai. 
M = 0. This implies that 



To compute A Q , we evaluate Eq. (p3|) at e 
lb = 1/6* from Eq. (|34l). In this case the integrals are all trigonometric 
integrals, and we have 

A = TT — aresin (b* l Q ) — aresin (b* l). 

Using Eq. (p9[) to express A in terms of ipobs, and making a small angle 
approximation, A is given by 
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A, 



VT - ?pobs ~ 



j>obs i 
lo ' 



The first Order term is given by 
A 1 = e 



~ d~ 








de 


e=0 V l ° 



b*Vi dl' r la b*Vi di' 



A 



A 



(35) 



(36) 



where 



A = ^2V2e b* 2 U 3 -b* 2 l l l2 + l. 



The derivative acts on both the e dependence in the integrals and in the 
upper limit £&, and care must be taken so that there is a cancellation of two 
divergent pieces which appear. Using a small angle expansion in ipobs, the 
first order correction to 6 is 



AV2M l 



obs 



Inserting the forms of A and Ai into Eq. (|3^) gives 

l 2V2R S l Q 

e = ifj obs (i + -) , 

vo Wöbs 



(37) 



(38) 



where Rs = 2M is the Schwarzschild radius. 

To lowest order, the physical distances in Figure [| are the inverse coor- 
dinate distances, 



l 



1 



1 



V2D LS V2D L 
and from Euclidean geometry, 9 is related to ß by 

ßD s 



e 



D 



LS 



Using these relationships in Eq. ( |38| ) and rearranging gives an approximate 
lens equation 



ß 



Dls + Di 
De 



obu 



2 R s Dls 
D s D L tp obs 
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which is the Standard result when D L + D LS = D s : 



01 2 ^5 D LS 1 

ß = ipobs =p-p: ; — • 39 

D s D L il) obs 

As a special case, we consider the Einstein Rings, an early prediction 
of "pre"-General Relativity only recently observed. If the source lies along 
the +z axis, directly opposite the lens from the observer at ß — © = 0, the 
observer sees the image as a circular ring surrounding the lens, or an Einstein 
Ring. In this special case, the lens equation, Eq. (|3Ö|) , is an implicit equation 
for the exact Observation angle for the ring: 



TT 



sin ip, 



obu 



lo y/JÖÖj Jl 
sin ijj obs 



dl' 



1 



dl' 



io 4m J h 



sin 2 ^ bs V* f(V) 
_ sin 2 V» bs l' 2 f(V) 

WWÖ) ) 



where l b is the postive root of Eq. (|3J) or the positive root of 
(sin tfj obs ) 2 l 2 b (l - 2V2Ml b ) -l 2 (l- 2V2MU) = 0. 



(40) 



(41) 



In terms of the future light cone of the source, an observer who sees the 
Einstein Ring is situated along the crossover line in Figure |l|. Points on this 
line are conjugate to the initial point, and the light cone has unstable singu- 
larities there. The crossover point in the cut function represents a limiting 
"Einstein Ring" at infinity, but the actual Observation angle for this ring is 
zero, so that the ring is not observable from infinity. 



5 Pseudo-Minkowski Coordinates 

As a final application of the cut function, we show that the so called pseudo- 
Minkowski coordinates form a well defined, global coordinate System for 
the Schwarzschild space-time with a constant density dust inferior. In this 
section, we use the complex stereographic angles (£, Q as coordinates on the 
sphere. 
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The pseudo-Minkowski coordinates are defined by integrals over the sphere 
at infinity of the cut function weighted against the first four Yj m , 

Xl , m = f Z(x a , C) y,, m (C) dS 2 = f l>m (xt), (l = 0, 1) (42) 
j s 2 

where 

2 dC A dC 



dS 2 



is the volume element on the sphere of the null generators of X + . There is a 
conceptual problem with the definition of the pseudo-Minkowski coordinates 
as stated in Eq. (|42|). Namely, the definition is ambiguous because the cut 
function u = Z[x a , £), is, in general, not Single valued at Z + , and so one does 
not know which portion of the cut to integrate over. 

The ambiguity is resolved by using the light cone structure to pull the 
integral back to the sphere of initial null directions at the initial point x a Q . 
To pull the integral back, we must have a function, 

( = C(<,V,V), (43) 

which relates the final angular positions at T + , the ((, () to the initial direc- 
tion of the geodesic, {rj,ff) at the initial point. Given a function of the form 
of Eq. (|4*3D, we can form the determinant of the Jacobian matrix, 

|J| = 9C9C_9CöC (u) 
di] dfj dfj drj ' 

and transform the integral from an integral over the sphere at null infinity 
into an integral over the sphere of initial directions: 



Xl , m = Z(x a ,C(x a ,v)) Y Wm ((K,v)) \ J\ dS 2 =f lm {<) = 0,1). 

J S Q 

(45) 

This integral defines the pseudo-Minkowski coordinates. 

We would like to show that the pseudo-Minkowski coordinates form a 
good coordinate System by showing that the Jacobian of the coordinate trans- 
formation defined by Eq. (|45"|) , 

Xlm = flm{x a o ), (46) 
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is non-zero. 

Due to the spherical symmetry of the space-time and the fact that the 
%i,m — ( x i,i-> x i,o-> ^1,-1) transforms as an 0(3) vector under space-time 
rotations, we can conclude that the functional form of the pseudo- Minkowski 
coordinates must be 

xi-i — X — iY — f(u , l , b, 7) sin 9 e"^, 
xi fi = Z = f(u , l , b, 7) cos 9, 
x hl = X + iY = f(u , h, b, 7) sin 9 e^, 

x 0fi = T = g(u ,l ,b,-f). (47) 

To test the non-vanishing of the Jacobian, all we need to do is to take a 
point of the (u , l ) plane, for example O = 0, and 9 = n, and check the 
transformation 

x ,o = g(u ,l ,b,i), x lfl = f(u ,l ,b,r/), (48) 

since this part of the coordinate transformation represents the "non-rotational" 
part. The determinant, D, of interest is given by 

= dxpo dx 10 _ dxpo dx w 
du Q dl dl du Q ' 

For points along the —z axis, using b and 7 as the initial parameters 
and the Jacobian expressing their relationship to the angles (C,C) found in 
Appendix B, the pseudo- Minkowski coordinates are 

f 39 

xi, m = J db A c?7 sin 6*00 (/ ,&)^^ u(u Q ,l ,b) F; jm (C(6, 7)), (1 = 0,1) 

(50) 

where the ränge in 7 is zero to 2tt and the ränge in b runs fully over both 
sheets of Solutions. The Integration over 7 does not cause any trouble for any 
of the integrals. At first glance, the convergence of the Integration over b is 
not clear, due to divergences in term d9 co /db as b approaches its maximum 
value, 
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h 1 



These divergences are all of order (b m — 6) _/3 with ß < 1, which ensures that 
the b integral also converges. The derivatives in question can be written as: 

dx 00 r f „ . a n ,v 00«, 

~dv~ = y sm Poo(io,o) 

^ ~dT J dbu ( l °' u °'fy s'm9 00 (l ,b) 



— -H = y rf sinö 00 (Z ,fe) cosöoo^o,?») 

— !- = V^37r — — / dbu(l ,u ,b) sin6' 00 (/ ,6) cos^^o, 6) — ^ 
3tt 9 /■ „ „ , x d 



9 i; J dbu(l ,u ,b) -((sin9 00 (l ,b)) 2 ). (51) 

The integrals are defined piecewise along the various segments of u(l Q , u Q , b) 
and 6*00 (/ , 6). The ränge in b runs from 5 = 0, when the null geodesic is ra- 
dially outgoing, and hence 9oo(l , b = 0) = ir, to a maximum value to 6 = b m , 
and, on the second sheet, back down to 6 = for radially ingoing rays, where 
9oo(lo, b = 0) = 0. The first integral is easily performed: 

'du' = ^ J db sinöoo ^°' & ) = ^ J d6 — (-cosöoo) 

= — v / tt(cosO — cos7r) = — 2y / 7r. (52) 

Likewise, 

dx W r- — f ■ n n i\ n n n ^oo 



. )it VSn j db sin (/ , 6) cos6» 00 (/ , 6) ; )J 

= ^ I db Tb{l {sin6c 



|^ ((sinO) 2 - (sinvr) 2 ) = 0. (53) 
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Therefore, when the initial point lies along the —z axis, the Jacobian of the 
transformation simplifies to 

D _ dx 00 dx w _ dxpo dxw _ _ 2 ^dxio 
du a dl Q dl du Q dl Q 

= -V3tt J db u(l ,u a ,b)-^ ((sinöoo^o,^) 2 ) , 

or finally 

D = -v^vrA X(/ ), (54) 

with 

l(l ) = J dbu(l ,u ,b)-^ ((sinöoo^o^)) 2 ) . (55) 

Thus, to determine if the pseudo- Minkowski coordinates are a good coor- 
dinate System, we only have to show that the integral T{1 ) has no extremum 
as the initial radial coordinate parameter, l a , is varied. An extensive numer- 
ical calculation shows that there are no extremum to this integral, whose 
values are plotted for many initial positions in Figure[5|. In fact, the integral 
is a constantly decreasing function, whose derivative is finite at all points l 
except l = 0, which corresponds to spatial infinity. Since spatial infinity is 
not a point in the space-time, we claim that the determinant, D, has a finite 
positive value for all initial positions. 

We have shown that the Jacobian of the transformation between the 
(u ,l ) and (xo,o,^i,o) portions of the total coordinate transformation, 

is non-zero. Using the spherical symmetry of the space-time and the inherent 
transformation properties of the x^ m , we can claim that the entire transfor- 
mation is non-singular, or that that the pseudo- Minkowski coordinates are a 
good coordinate System of Schwarzschild space-time with a constant density 
dust inferior. 
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A The Initial Direction and the Parameter b 



The parameter b, which arose as a constant of Integration when integrating 
the null geodesic equations, parametrized the initial direction of the geodesic. 
In this appendix, we choose the motion of the geodesic to remain in the x-z 
plane and the initial point to lie on the —z axis. The initial direction of a 
geodesic is captured by giving an angle, ip, between the spatial part of the 
directed tangent vector to the geodesic and the z axis. We are interested in 
determining the relationship between the angle 1p and the parameter b. 

From Eq. (|TT|) , the coordinates of a null geodesic restricted to the x-z 
plane were given in terms of / by 



u = u + 



1 ± yjkjl') A(l>) 

A{V){21>* f{l')) j 




(56) 

with 

A(V\ = - 

f(l>) h{V) 

Up to rescaling, the (null) tangent vector to this geodesic is 

(du <M de d£\ (du dß_ \ 
\dVdVdVdl) \dV 'dV )' [ ' 

Since both the measure of angles and the length of a null vector are 
independent of conformal factors, any conformally related metric may be 
used to compute them. In what follows we use the physical metric of our 
model, which in the coordinates (t, l = l/(V2r), 9, <fi) is 

ds 2 = f{l)dt 2 - ^dl 2 - ^dü 2 = f{l)dt 2 - gij dx l dx J , (58) 

where /(/) and h{l) are the coefficients of the metric given in Eq. (||) and 
is a spatial metric. The spatial part of null vector, Eq. (|57|), normalized in 
the physical metric, is the three vector 
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with 



L 



\ 2l 4 + 2/ 2 \ dl ) 



h(l) 1 (d6\ 2 



The value of the derivative d9/dl is determined using Eq. (|56|): 

(d9\ 2 _ f{l)h{l)b 2 

\di) i-b 2 i 2 f{iy 



(60) 



A unit spatial vector pointing in the radial direction is given by 



( 



s/W) 



V21 2 



0,0 



The inner product between L l and r % gives the angular direction of the 
geodesic, namely, 



After some algebra, Eq. flBTD can be solved for b, giving our desired result 



The ränge in b at the initial point is determined by Eq. (f32|) . The param- 
eter ranges from b = 0, corresponding to radially outgoing rays when 1p = 0, 
to a maximum value, b = b m , where iß = n/2, back down to b = where 
ip = Ii, and again the geodesic travels radially. Either b or iß may be used to 
parametrize the initial direction of the geodesic. 

B Füll Angular Dependence of the Light Cone 

In Section 2, we integrated the null geodesics emanating from a point on the 
— z axis, restricted to the x-z plane, in terms of a parameter l. The angular 
integrals were 



gijf l L J = cos iß. 



(61) 



b 



(62) 
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= 0, 




Q = e(l,lo,b)=ir- dl' ±^= V (63) 



We want to perform a rigid rotation of this restricted Solution restoring the 
füll angular dependence, and allowing the initial point to be at any position. 

Due to spherical symmetry, the geodesic equations separate into two 
time/radial equations and two angular equations. An arbitrary Solution to 
the angular part of the geodesic equations can be obtained by performing a 
rigid rotation of the Solution given in Eq. (|63|). For such a Solution, the mo- 
tion will take place in a new plane, but the angle 0(/, Z D , b) will be preserved. 

To perform the rotation we use complex stereographic coordinates, ((, (), 
as coordinates on the sphere defined in Eq. (|T2|) . In terms of C, the Solution 
corresponding to Eq. (Iii]) is 



C(J, lo, b) = cot e( ^°' b) = C(J, lo, b). (64) 



Under an SU{2) rotation, Q transforms as 



a C + b acot^ + b 
C < + d ccot^ + d' [ ' 



where a, b, c, d are the Cayley-Klein parameters |p.Q|| , which can be expressed 
in terms of Euler angles, a, ß, and 7 as 



a = sm — e 2 
2 



b = cos-e^-^. 



2 

c = — cos — e2 a p > 
2 

d = mn|e^-T-^. (66) 

To determine the values of the Euler angles, we note that when 6 = tt, the 
geodesic is at the initial position, ( Q . From C' — Co; we have: 
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^ = cot -e i/3 = cot (67) 
d 2 2 

This condition fixes two of the Euler angles, a and ß, to be a = 9 a and 

ß = 4> - Thus, in terms of the new initial point, Co, the Cayley-Klein 

Parameters are 



a 



b 



c 



d 

The remaining free parameter 7 gives the orientation of the plane in which 
the geodesic moves. In the case that the initial point lies on the — z axis, 7 is 
the angle 0. When the initial point of the geodesic is rotated to an arbitrary 
location, the parameter 7 acts as an angle about the new axis of symmetry 
in the System. 

Our final, füll Solution to the angular part of the geodesic equations is 
obtained using Eqs. (|68[) with Eq. (|6o|): 



1 



'Co 

1 + C0C0 vco, 



C0C0 



'Co 



\ 1 + C0C0 VC 



Co Co I Co 



\ 1 + C0C0 VCo, 



1 



'CoV 



1 + C0C0 VCo, 



— e 2 



(68) 



V C °/ V-VCZe i7 cot^^ + e-^; 

where we have dropped the prime on C- The angular Solution, C is a function 
of the initial point (ZojC^Co)? a parameter along the light cone l, and two 
free parameters, (b, 7), which span the sphere of initial null directions at the 
initial point. The dependence on l , l, and b comes through Q(l,l ,b) by 
integral expression 



/ ■ 1 b <fk(jj , 

e(z, i , b) = TT - / di' ± - v (70) 
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When the value of / is taken to zero, Eq. (|69|) gives the final angular 
location of a point on I + in terms of initial directions (b, 7) and the initial 
point x" = (u , l , Co, Co)- In this case, we denote C(^o, Co, Co, 1 = 0, b, 7) by 
Coo(^o, b, 7), and Q(l = 0, l , b) by Öqq. The existence of such a function 
provides a mapping from the sphere of initial null directions, (b, 7), to the 
sphere of null generators at X + . The Jacobian matrix of the mapping is given 



In Section 5, we use the determinant of this Jacobian to transform integrals 
over X + to integrals over the initial null directions. 
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Figure 1: The light cone of a point which has formed conjugate points in its 
future. The crossover line represents points in the space-time where Einstein 
Rings are observed, while the line of caustics are stable singularities conjugate 
to the initial point. 
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Figure 2: Three cut functions with one dimension suppressed (5 x R). As 
the initial radial position moves away from the spatial center smooth cuts 
give way to cuts with singularities. 
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Figure 3: A cusp ridge appears on the light cone cut from an initial point 
sufficiently far away from the center of the space-time. Near X + , the intersec- 
tion of the light cone of this initial point with a time-like surface would look 
similar to this cut, except that the "umbrella" at the top of the cut would be 
inside the main body of the wavefront. 
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Source 




Observer 

Figure 4: The schematic representation of the path of a geodesic observed 
in gravitational lensing. Distances between the lens and observer, lens and 
source, and observer and source are shown along with the Observation angle, 
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Figure 5: A plot of the integral, X(/ Q ), as a function of the initial radial 
Position shows that there are no extrema. The integral has a finite derivative 
for all points except l = 0, which is not in the physical space-time. 
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